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INTRODUCTION 


The requirements placed on an unsteady aerodynamic theory intended for turbomachinery aero- 
elastic or aeroacoustic applications will be discussed along with a brief description of the various 
theoretical models that are available to address these requirements. The major emphasis is placed on 
the description of a linearized inviscid theory which fully accounts for the affects of a nonuniform 
mean or steady flow on unsteady aerodynamic response. Although this linearization has been 
developed primarily for blade flutter prediction, more general equations will be presented which 
account for unsteady excitations due to incident external aerodynamic disturbances as well as those 
due to prescribed blade motions: In this presentation we will focus on the motivation for this 

linearized unsteady aerodynamic theory, outline its physical and mathematical formulation and 
present examples to illustrate the status of numerical solution procedures and several effects of 
mean-flow nonuniformity on unsteady aerodynamic response. This presentation is based on a paper 
of the same title which is published in full in the Proceedings of the Tenth U. S. National Congress 
of Applied Mechanics (ref 1). 


• Linearized unsteady aerodynamic analysis 

• Real blade geometry 

• Mean blade loading 

• Shock phenomena 

• Design applications 

• Aeroelastic 

e Blade flutter and forced vibration 

• Aeroacoustic 

• Noise generation, transmission 
and reflection 
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MAJOR ASSUMPTIONS 


The development of theoretical models to predict unsteady flows through turbomachines is a 
formidable task. The analyst is confronted with determining the time-dependent, three-dimensional 
flow of a viscous compressible fluid through a geometric configuration of enormous complexity. This 
task has required the introduction of a considerable number of simplifying assumptions to make 
the problem mathematically tractable and to render the resulting solutions useful to designers. 
For the most part, the theoretical formulations that have been developed to predict the unsteady 
aerodynamic phenomena associated with blade flutter or forced vibration consider the blades of 
an isolated two-dimensional cascade, neglect viscous effects at the outset and regard unsteady 
fluctuations to be of sufficiently small amplitude so that a linearized treatment of the unsteady 
perturbation is justified. In addition the resulting two-dimensional inviscid flow is assumed to 
remain attached to the blade surfaces, the mean flow is assumed to be at most a small irrotational 
steady perturbation from a uniform stream at the cascade inlet, and any shocks that might occur 
are assumed to be of weak to moderate strength and have small curvature. 


• Isolated blade row 

• Two-dimensional inviscid (Re-^oo ) flow 

• Small-amplitude periodic unsteady 
excitations 

• Attached flow 

• Irrotational mean flow at Inlet: 

7= v m00 + v? 

• Weak shocks 
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REQUIREMENTS 


In general, the unsteady aerodynamic analyses intended for turbomachinery aeroelastic appli- 
cations must be applicable to fan, compressor and turbine cascades, to subsonic, transonic and 
supersonic Mach numbers and to moderate through high frequency structural and external aero- 
dynamic excitations. Then, to determine the aeroelastic and aeroacoustic characteristics of the 
blading such analyses must be capable of predicting the unsteady loads acting on the blades and 
the amplitude and wave numbers of the acoustic waves which carry energy away from the blade 
row and the entropic and vortical fluctuations which are convected downstream. These responses 
arise from the various sources of unsteady excitation including prescribed blade motions, variations 
in total temperature and pressure (“entropy and vorticity waves”) at inlet, and variations in static 
pressure (acoustic waves) at inlet and exit. For blade flutter applications it is only necessary to 
predict the unsteady loads acting on the blades as a result of prescribed blade motions; for forced 
response applications the unsteady blade loads due to incident entropic, vortical and acoustic dis- 
turbances are also required. Finally, for aeroacoustic applications the parameters associated with 
far-field acoustic responses must be determined. 


• Fan, compressor and turbine cascades 

• Subsonic, transonic, supersonic Mach numbers 

• Moderate to high excitation frequencies 

• Response predictions 

• On blades: surface pressures, global unsteady airloads 

• Far field: outward propagating acoustic waves 
vorticity and entropy variations downstream 

• Prescribed excitations 

• Blade motions (flutter) 

• External aerodynamic disturbances (forced vibration) 

• Vortical and entropic disturbances at inlet 

• Acoustic disturbances at inlet and exit 


i 
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TRANSONIC CASCADE: < M_ 0 0 < 1 


A representative cascade configuration is shown in the figure below which depicts a two-dimen- 
sional section of a transonic compressor rotor ( M < M_ 00 < 1. The cascade stagger angle is 
denoted by 0 and the blade spacing by G. In the absence of unsteady excitation the blades are 
identical in shape, equally spaced and their chord lines are oriented at the same angle, 0, relative 
to the axial flow direction. The inlet and exit free-stream flows are described by the velocity vectors 
V^poo. The free-stream flow angles measured relative to the axial-flow (or £-) direction are denoted 
by For the configuration illustrated the inlet and exit conditions are such that normal shocks 

(Sh) emanate from the blade suction surfaces and vortex wakes (W) emanate from the blade trailing 
edges and extend downstream. 
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UNSTEADY EXCITATIONS 


The unsteady fluctuations in the flow arise from one or more of the following prescribed sources: 
blade motions, upstream and/or downstream acoustic disturbances which carry energy toward the 
blade row, and upstream entropic and vortical disturbances which are convected through the blade 
row by the mean flow. These excitations are assumed to be of small amplitude and periodic in time. 
The external aerodynamic excitations are also spatially periodic, while the structural excitation 
is periodic in the “circumferential” or 77-direction. For example, we consider blade motions and 
incident acoustic disturbances as described below. Here 7 ? measures the displacement of a point on 
a moving blade surface relative to its mean or steady-state position, X is a position vector, m is a 
blade number index, t is time, is the reference-blade (m = 0) complex displacement-amplitude 
vector, u> and <7 are the temporal frequency and interblade phase angle, respectively, of the unsteady 
excitation and Re{ } denotes the real part of { }. Also, p/ )TO o and ic Too are the amplitude and wave 
number, respectively, of an incident pressure fluctuation, pj(X,t), coming from far upstream (—00) 
or far downstream (+00). Note that the interblade phase angle, <7, of an incident disturbance is 
The temporal frequency and wave number of an incident vortical or entropic disturbance 
are related by = — « Too • G, but a more complicated relationship exists between u> and « ?0O for 
an incident pressure disturbance. 



Blade motions: R (X + mG, t) = Re j r*(Xj exp [i (cot + mo)] J 
Incident disturbances: p'- «> (X, t) = Re jp^ooexp [i(/c_ oo* X + cot)] | 
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TIME-DEPENDENT GOVERNING EQUATIONS 


The equations governing the fluid motion follow from the integral forms of the mass, momentum 
and energy conservation laws and the thermodynamic relations for a perfect gas. The former provide 
a coupled set of corresponding nonlinear differential equations (the Euler equations) in continuous 
regions of the flow and jump conditions at surfaces across which the inviscid flow variables are 
discontinuous, i.e., at vortex-sheet wakes and shocks. In continuous regions the energy equation can 
be replaced by the requirement that the entropy following a fluid particle must remain constant. 
In addition to the foregoing field equations and jump conditions, the attached flow assumption 
requires that the unsteady flow must be tangential to the moving blade surfaces and information 
on the uniform inlet and exit flow conditions and the incident entropic, vortical and acoustic or 
static pressure disturbances must be specified. The remaining steady and unsteady departures from 
the uniform inlet and exit conditions must be determined as part of the time-dependent solution. 
This foregoing aerodynamic problem is a formidable one as it involves a system of nonlinear time- 
dependent equations with conditions imposed on moving blade, wake and shock surfaces in which the 
instantaneous positions of the wakes and shocks must be determined as part of the solution. Because 
of these features and the prohibitive expense that would be involved in obtaining the aerodynamic 
response information needed for aeroelastic or aeroacoustic applications, the usual approach is to 
examine limiting forms of the full governing equations with the intention of providing efficient 
analyses for design applications. 


• Integral conservation laws 

• Thermodynamic relations 



• Euler equations at field points 

rsj 

• Jump conditions at moving shocks (Mf ± 0) and at 

vortex sheet boundary layers and wakes (Mf = 0) 

• Flow tangency condition at moving blade surfaces 



• Far-field behavior 
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UNSTEADY AERODYNAMIC LINEARIZATIONS 


Because of the complexity of the nonlinear time-dependent unsteady aerodynamic problem, 
linearized treatments of the unsteady flow are often considered. The major linearizations that 
have been proposed are the following: classical linearized theory, time-linearized transonic flow 
theory and the present theory in which unsteady disturbances are regarded as small relative to 
a fully nonuniform mean flow. The essential differences between these theories arise from the 
manner in which the steady flow is represented. In classical theory both steady and unsteady 
departures from a uniform stream are regarded as small and of the same order of magnitude. 
In time-linearized transonic theory steady and unsteady disturbances are regarded as small and 
very small, respectively, relative to uniform free-stream flow properties. Finally, in the present 
linearization no restriction is placed on the steady flow but the unsteady perturbations are assumed 
to be of small amplitude. The classical theory applies at the (reduced) frequencies of interest for 
turbomachinery applications, but steady flow variations have no impact on the unsteady response. 
Time-linearized transonic theory applies at Mach numbers near one and the unsteady perturbation 
depends on the steady flow, but this theory is formally restricted to low-frequency unsteady motions. 
The present theory fully includes the effects of nonuniform mean flow and applies throughout the 
Mach number and frequency range of interest for turbomachinery applications. 


• Classical theory 

P U,y,t) = P-oo + p(x,y) + Re | p(x,y)e lwt | + •• 

• Time-linearized transonic theory 

p ( *i y «t ) * P-® + P ( x,y ) + Re |p(x,y)e ia ' t i+... 

• Present theory 

P (x,y,t) = P(x,y) + Re j p(x,y)e lwt ( + . . . 
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PRESENT LINEARIZATION 


The equations governing small-amplitude unsteady departures from a nonuniform mean flow are 
determined by expanding the various flow variables in asymptotic series in e, where e is a measure of 
the amplitude of the unsteady excitation. Thus, for example, in the first equation below P(X,t) is 
the time-dependent fluid pressure, P(X) is the pressure in the steady background flow, Re{p(X)e lujt ] 
is the first-order time-dependent pressure and p(X) is its complex amplitude. In addition, Taylor 
series expansions and surface vector relations are used to refer information at a moving blade, wake 
or shock surface (<S) to the mean position of this surface (S). Equations governing the zeroth-order 
or steady and the complex amplitudes of the first-order unsteady flow properties are obtained after 
substituting the foregoing expansions into the full time-dependent governing equations, equating 
terms of like power in e and neglecting terms of higher than first order in e. It follows from the 
original assumptions that the steady background flow is governed by a full-potential boundary- 
value problem and that the complex amplitudes of the unsteady flow properties are governed by 
a system of time-independent linearized equations with variable coefficients which depend on the 
underlying mean flow. In the unsteady problem surface conditions can be imposed at the mean 
surface locations, and in both the steady and first-order unsteady problems the required solution 
domain can be limited to a single extended blade- passage region. 


• Series expansions 

P (X,t)= P(X) + Re {p(X)e ia,t }+ ... 



(P + /? .yP + ...) s 



fg* f (t s ,/? ) + 


• Full potential boundary-value problem for steady 
background flow 

• Linear variable-coefficient boundary-value problem for first- 
order unsteady flow 

• Time-independent 

• Surface conditions at mean surface locations 

• Single extended blade-passage solution domain 


I 
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THE STEADY BACKGROUND FLOW 


As a consequence of our assumptions regarding shocks and the steady flow far upstream of 
the blade row, the mean or steady background flow through the cascade will be homentropic and 
irrotational; i.e., V = where V and $ are the local steady velocity and velocity potential, 
respectively. The field equations governing the steady flow follow from the mass and momentum 
conservation laws and the isentropic relations for a perfect gas. Here, p, P , M and A and are the 
local steady density, pressure, Mach number and speed of sound propagation, respectively, and 7 is 
the specific heat ratio of the fluid. Surface conditions for the zeroth-order or steady flow apply at 
the mean positions of the blade ( B ), wake (W) and Shock (Sh) surfaces. Blade mean positions are 
prescribed, but the mean wake, i.e., the stagnation streamlines downstream of the blade row, and 
shock positions must be determined as part of the steady solution. Since, by assumption, the flow 
remains attached to the blades, a flow tangency condition applies at blade surfaces. In addition 
the steady pressure and normal velocity component must be continuous across blade wakes and 
mass and tangential momentum must be conserved at shocks. Finally, three of the far-field uniform 
velocity components, or the equivalent information, must be prescribed to completely specify the 
steady boundary-value problem. The fourth or remaining component can be determined in terms 
of the three prescribed using the integral form of the mass conservation law. 


• Field equations : V(/>VS>) s O 

p {r ~' ] ’(rMljP/2) (r ' ll/y =(M. 00 A) 2 *F|V<H 

• Surface conditions 

• Blades (B): V<J>1T*0 

• Wakes (W):[V$]-n s O,ll<I>J = constant 

• Shocks (Sh):8>V<t«]| n* |4>1*0 

• Far-field conditions 

• Uniform flow conditions 

• Analytic solutions for steady disturbances 
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THE LINEARIZED UNSTEADY FLOW - 1 


The field equations governing the first-order unsteady perturbation of a nonlinear homentropic 
and irrotational steady flow are determined from the full time-dependent mass, momentum and 
entropy transport equations and the thermodynamic equation relating the entropy, pressure and 
density of a perfect gas. These equations can be cast in a convenient form by introducing the 
Goldstein (ref. 2) velocity decomposition, i.e., V = v R + S/<f>. The rotational component of the 
unsteady velocity (u H ) is divergence-free far upstream of the blade row and it is independent of 
the pressure fluctuation (p); the irrotational component (S/(f>) * s related directly to the unsteady 
pressure fluctuation. The resulting field equations for the first-order entropy (.s), rotational velocity 
(v R ) and velocity potential (cf>) are given below. First-order partial differential equations describe 
the transport of entropy and rotational velocity through the blade row. The unsteady potential 
is governed by a second-order equation, which is locally elliptic at field points at which the local 
steady Mach number (M) is less than one and locally hyperbolic at those points at which M > 1. 
Note that the rotational velocity provides a forcing function term to the potential equation. Also, 
if there are no entropy and rotational velocity fluctuations at inlet, then only a single field equation 
must be solved to determine the unsteady flow. 

• Velocity decomposition .• 7 = 7 R ♦ v<#> 

where T R = f (s-^v^VS)) and p=-p ^ 


• Field equations 



^-(v r -sv<!>/ 2) ♦ [(v R - sv<J>/2)-v] v4> a 0 
Dt (A " 2 of >-?"'*• •£v R) 
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THE LINEARIZED UNSTEADY FLOW - 2 


Conditions on the linearized unsteady perturbation at Blade ( B ), wake (IT) and shock (Sh) 
mean positions are obtained in a similar fashion, i.e., by substituting the asymptotic, Taylor and 
surface-vector series expansions into the full time-dependent surface conditions and equating terms 
of like order in e. The resulting first-order flow tangency, wake-jump (continuity of normal velocity 
and pressure) and shock-jump (conservation of mass and tangential momentum) conditions are 
indicated schematically below. Note that the blade displacement (rg) is prescribed but that the 
normal(to the shock) component of the shock displacement (r*sh • n) must be determined as part 
of the unsteady solution. Wake displacements have no 'impact on the linearized unsteady problem. 
In addition to the surface conditions, we require information on the unsteady flows far upstream 
and far downstream from the blade row. In these regions the linearized unsteady equations reduce 
to the constant coefficient equations of classical linearized theory, and analytic solutions for the 
velocity potential fluctuations (cf> p ) due to acoustic response disturbances and the far-downstream 
potential fluctuations (<f> R ) associated with the vortical or rotational velocity disturbances can be 
determined. These analytic far-field solutions can be matched to near-field numerical solutions, and 
they thereby serve to complete the specification of the linearized unsteady boundary-value problem 
(ref. 3). 


• Surface conditions 

• Blades: <f>„ * f (4>,7 B ,v R ) 

•Wakes: ] * [pJ =0 


• Shocks: |[F(<J>)«#> + G (*)<#>„] = f (*,7^ n .v R ) 

•sh,n * ” 1^1/ IM, I s ]l * [v R J =0 

• Far -field conditions 

s -co »^-a> Prescribed 


* 




*,N 


+ +0UT + * 


Prescribed 
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NUMERICAL SOLUTION DOMAIN 


The foregoing steady and linearized unsteady boundary-value problems account for the effects 
of blade geometry, mean blade loading and transonic, including moving shock, phenomena on the 
unsteady fluctuations arising from small-amplitude harmonic excitations. The unsteady equations 
are linear, time-independent and contain variable coefficients which depend on a fully nonlinear ho- 
mentropic and irrotational steady background flow. Numerical resolutions of the nonlinear steady 
and the linearized unsteady problems are required to determine the aerodynamic response informa- 
tion needed for aeroelastic and aeroacoustic applications. Because of the cascade geometry and the 
assumed form of the unsteady excitations (i.e., periodic in t and 77), such resolutions are required 
only over a single extended blade-passage region. In addition, since analytic far-field unsteady so- 
lutions have been determined, the numerical solution domain can be restricted further to a single 
extended blade-passage region of finite extent in the axial direction as shown below. 




W: Continuity of pressure and normal 


velocity across wakes 


Sh: Conservation of mass and tangential 


momentum across shocks 


P: Periodicity upstream 
M: Near/far field matching 
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AERODYNAMIC RESPONSE AT A BLADE SURFACE 


For aeroelastic and aeroacoustic applications, solutions to the nonlinear steady and the linearized 
unsteady boundary-value problems are required to predict the unsteady aerodynamic response at 
a moving blade surface (i.e., the unsteady surface pressures and the global unsteady airloads) and 
in the far-field (i.e., the unsteady pressure fluctuations), respectively. In particular, the pressure 
acting at the instantaneous position (5) of a given blade surface is made up of two components: a 
harmonic component (px) which is determined by the steady ($) and the linearized unsteady (<f>) 
potentials and the prescribed blade displacement (fa), and an anharmonic component (p Sh ) which 
is caused by the motion (rsh ,b) of a shock along the blade surface. The anharmonic pressure is 
determined by analytically continuing the solution to the steady boundary- value problem from the 
mean to the instantaneous shock location (ref. 4). Although the pressure disturbance ps is not 
everywhere harmonic, its regions of an harmonicity are small. Consequently, the first order global 
unsteady airloads are harmonic in time (ref. 5). In particular, if each two-dimensional blade section 
undergoes a pitching oscillation (re = a x R v ) about an axis fixed to the blade, the first-harmonic 
unsteady moment is determined by integrating the product of the first-harmonic component of the 
unsteady surface pressure and R p • f over the mean blade surface and subtracting a term consisting 
of the product of the steady pressure jump across the shock, R v • r and the shock displacement 
along the blade surface. 


Surface pressure 

P^= P H ($, 0. n> + P Sh <$- r Sh, B’tH - 



Unsteady moment 

m = /PHRp* d ^- r sh,B I t P B^P ,ir )lB + - 
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NUMERICAL RESULTS 


At this point we have completed our description of the unsteady aerodynamic model and pro- 
ceed to present numerical results to partially illustrate the status of numerical procedures for solving 
the nonlinear steady and linearized unsteady boundary-value problems and to demonstrate several 
important effects associated with nonuniform steady flow on the aerodynamic response at a moving 
blade surface. We refer the reader to refs. (6 and 7) for a description of the numerical procedures 
used. We will present results for two-dimensional compressor- or fan-type cascades operating at 
subsonic inlet and exit conditions. Theoretical results for steady surface Mach number (M) distri- 
butions, first-harmonic unsteady pressure-difference (A ph) distributions and unsteady aerodynamic 
moments (m) will be presented for blades undergoing pure pitching (torsional) motions with a = 1 , 0 
about their midchords. The stability of such motions depends upon the sign of the out-of-phase 
moment (m/). If m/ > 0, the airstream supplies energy to the blade motion, and this motion is 
unstable according to linearized theory. We will consider a subsonic cascade of NACA 0012 air- 
foils to illustrate the effects of a relatively thick, blunt-nosed blade geometry and variable mean 
incidence on the unsteady response and a subsonic/transonic cascade of 5% thick flat-bottomed 
double-circular-arc (DCA) airfoils to illustrate the effects of mean blade loading and transonic phe- 
nomena on the response at high subsonic inlet Mach number. For purposes of comparison results 
for flat-plate cascades, operating in uniform mean flows will be included along with those for the 
NACA 0012 and DCA cascades. The example cascades each have a stagger angle 0 of 45 deg and 
a unit gap/chord ratio ( G = 1). The steady flows through the NACA 0012 and DCA cascades have 
been determined by imposing a zero-load condition at sharp blade edges. 


• Torsional vibrations about midchord: r B =(ce x Rp) 

• Response parameters: Ap H and m 

• Subsonic flow: NACA 0012 airfoils 

• Subsonic/transonic flow: DCA airfoils 

mu** 
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SUBSONIC FLOW: EFFECT OF INCIDENCE - 1 


Steady and unsteady flows through the staggered cascade of modified NACA 0012 airfoils have 
been determined for an inlet Mach number (M_ oc ) of 0.6 and four inlet flow angles. The predicted 
surface Mach number distributions for fl-oo = 48, 50, 52 and 54 deg are shown below. The calculated 
exit Mach numbers are respectively 0.595, 0.557, 0.522 and 0.490, and in each case the calculated 
exit flow angle is approximately 47.7 deg. These steady flows are entirely subsonic with a peak Mach 
number of 0.789 occurring at x = 0.113 on the pressure (lower) surface of the blade for = 48 
deg, and 0.8, 0.86 and 0.96 occurring at x = 0.07, 0.05 and 0.03 on the suction (upper) surface for 
Cl-oo = 50, 52 and 54 deg, respectively. In each case the mean flow stagnates within 0.2% of blade 
chord downstream from the leading edge. 
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SUBSONIC FLOW: EFFECT OF INCIDENCE - 2 


Unsteady response predictions for the NACA 0012 cascade and for a flat-plate cascade operating 
at = 0.6 and fl _ <*, = 45° are shown on this and the next twofigures. Shown beloware unsteady 

pressure-difference distributions and aerodynamic moments for the reference (m = 0) NACA 0012 
and flat-plate blades undergoing unit-frequency pitching motions at a = 90 deg. The unsteady 
pressure difference is singular and behaves like a multiple of x ~ near the leading edge of the flat- 
plate airfoil. In contrast, the unsteady pressure is analytic in the vicinity of the rounded leading 
edge of the NACA 0012 blade. In this case both the real and imaginary components of the unsteady 
pressure difference are zero at the leading edge and reach local extrema very close to the leading 
edge. The results indicate that the coupling between the steady and unsteady flows, due to blade 
geometry and mean loading, leads to a reduction in the out-of-phase pressure difference, Im{Ap(x)}, 
over a forward part of the NACA 0012 blade and, therefore, a reduction in the out-of-phase moment 
opposing the blade motion. 


Unsteady pressure difference distributions: NACA 0012 cascade; 

0 = 450, G = 1.0, M _ 00 = 0.6, co = 1.0, o = 90 °. 

NACA 0012 cascade, flat-plate cascade 
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SUBSONIC FLOW: EFFECT OF INCIDENCE - 3 


Unsteady moments acting on the reference blade of the NACA 0012 cascade operating at inlet 
flow angles of 48, 50, 52 and 54 deg and on the reference blade of the flat-plate cascade operating 
at Q-oo = 45 deg are shown below for blades undergoing unit-frequency torsional vibrations over 
the entire range of interblade phase angles, i.e., cre[— 7r,7r]. The abrupt changes in the moment 
curves are indicative of an acoustic resonance. The blade motions are superresonant (i.e., acoustic 
response disturbances persist in the far field and carry energy away from the blade row) at interblade 
phase angles lying between the lowest and highest resonant phase angles and subresonant (acoustic 
response disturbances attenuate with increasing distance from the blade row) at the interblade phase 
angles below the lowest and above the highest resonant phase angles. The blade motions considered 
below are stable but the NACA 0012 results indicate that the effect of mean blade loading tends 
to be destabilizing. Note that for a given a the out-of-phase moment moves closer to the stability 
boundary as the inlet flow angle is increased. 


Unsteady moment vs. interblade phase angle: 0 = 45°, G = 1.0, M_oo = 0.6. 
NACA 0012 cascade, flat-plate cascade 
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SUBSONIC FLOW: EFFECT OF FREQUENCY 


The effect of frequency on the out-of-phase component of the unsteady moment due to torsion 
about midchord is illustrated below for the NACA 0012 cascade operating at oo = 54 deg and 
for the flat-plate cascade operating at ft_oo = 45 deg. The NACA 0012 blades experience a region 
of subresonant torsional instability for u? = 0.25, 0.5 and 0.75, with the extent of this region 
decreasing with increasing frequency. The subresonant torsional motions of the flat-plate cascade 
are unstable only at the lowest frequency considered, i.e., u = 0.25. Thus the nonuniform flow 
through the NACA 0012 cascade extends the frequency range over which the blades are susceptible 
to a torsional instability. 


Unsteady moment vs. interblade phase angle: 
0 = 45°, G = 1.0, M _ oo = 0.6. 


NACA 0012 cascade, Q_oo = 54° Flat-plate cascade, Q_oo = 45° 
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SUBSONIC/TRANSONIC FLOW: EFFECT OF MACH NUMBER - 1 


We now consider the staggered cascade of sharp-edged double-circular-arc airfoils. In particu- 
lar, the airfoils have flat lower surfaces, circular arc upper surfaces and maximum thicknesses at 
midchord of 0.05. Full-potential steady and linearized unsteady flows through this example config- 
uration have been determined for inlet Mach numbers of 0.5, 0.7, 0.8 and 0.9. The steady flows 
have been determined by imposing a zero-load requirement at blade leading and trailing edges. As 
a consequence, only the inlet Mach number is prescribed with the remaining inlet and exit param- 
eters determined as part of the steady flow solution. Numerical results for this configuration along 
with those for a corresponding flat-plate cascade (0 = 45 deg, G = 1) are given in this and in 
the following three figures. Shown below are the predicted surface Mach number distributions for 
the example DCA cascade. For the prescribed inlet Mach numbers stated above, the calculated 
exit Mach numbers are 0.43, 0.57, 0.62 and 0.64, respectively. In addition, the calculated inlet flow 
angles are 49.0, 49.2, 49.4, and 49.6 deg, respectively, and in each case the calculated exit flow angle 
is approximately 43.0 deg. The steady flows at A/_oo = 0.5, A/_ <» = 0.7 and A/_oo = 0.8 are entirely 
subsonic with the maximum suction-surface Mach numbers of 0.561, 0.804 and 0.941 occurring at, 
respectively, 40.8, 38.5 and 36.5% of blade chord downstream from the leading edge. The steady 
flow at A/_ oc = 0.9 is transonic with the supersonic region extending from 18.5 to 52.5% of blade 
chord along the suction surface and terminating at a shock discontinuity. The Mach numbers at 
the foot of the shock are 1.193 on the upstream or supersonic side and 0.871 on the downstream or 
subsonic side. 


Surface Mach number distributions: DCA cascade; 0 = 45°, G = 1.0 
Suction surface Pressure surface 
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SUBSONIC/TRANSONIC FLOW: EFFECT OF MACH NUMBER - 2 

The effect of Mach number on the response to in-phase (<r = 0 deg) unit-frequency torsional blade 
vibrations is illustrated below for the DCA and flat-plate cascades. The blade motions at a — Odeg 
are superresonant. For the flows at M - = 0.5, 0.7 and 0.8 two acoustic waves persist in the far 
field — one upstream and one downstream— and propagate away from the blade row. For the DCA 
cascade operating at M_ <» = 0.9 there are three such waves — two upstream and one downstream. 
Finally, for the flat-plate cascade operating at M_ = 0.9 there are four such waves — two upstream 
and two downstream. The out-of-phase pressure-difference distributions and unsteady moments for 
the reference DCA and flat-plate blades reflect this change in character of the acoustic response in 
the far field, since the trends indicated by the results for M_oo = 0.5, 0.7 and 0.8 are not sustained 
at M-oo = 0.9. Also, a comparison of the DCA and flat-plate pressure-difference curves for in-phase 
motions suggests that the influence of mean flow gradients on the unsteady aerodynamic response 
becomes more pronounced with increasing Mach number. The pressure difference distributions 
for the DCA and flat-plate blades are very similar for the two lower inlet Mach number, differ 
somewhat for M _<*, = 0.8 and differ substantially for M_ «*, = 0.9. The differences at M_ = 0.8 
can be attributed to the relatively large gradients in the subsonic mean flow that occur along the 
suction surface of each DCA blade. The substantial differences at M_oo = 0.9 are caused by the 
transonic effects associated with the DCA cascade and by the different far-field acoustic response 
environments produced by the two cascades. 


Unsteady pressure-difference distributions and unsteady moments: 

0 = 45°, G = 1.0, CO = 1.0, 0 = 0°. 
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SUBSONIC/TRANSONIC FLOW: EFFECT OF FREQUENCY 


Shown below are first-harmonic, out-of-phase, pressure-difference distributions and aerodynamic 
moments for the example DCA and flat-plate cascades operating at an inlet Mach number of 0.9. 
Here the blades are undergoing out-of-phase (<r = 180 deg) torsional vibrations about midchord 
at different prescribed frequencies. Recall that for a discontinuous transonic flow there are two 
contributions to the first-harmonic unsteady moment: one arising from the harmonic unsteady 
surface-pressure response and the other from the anharmonic surface pressures produced by shock 
motion. However, for the DCA cascade at M_ 0 0 = 0.9, the mean shock location is only slightly aft 
of blade midchord and, therefore, the anharmonic surface pressures make only a small contribution 
to the unsteady moment. A comparison of the DCA and flat-plate results depicted below indicates 
the dramatic impact of transonic mean-flow phenomena on unsteady aerodynamic response. A 
second interesting feature indicated by these results is the change in the unsteady moment behavior 
as the blade vibration frequency is increased from 0.25 to 0.5. This change occurs because the 
out-of-phase blade motions of the DCA and flat-plate cascades are subresonant for u> = 0.1 and 
0.25 and superresonant for u> = 0.5, 0.75 and 1.0, and this change in the far-field acoustic response 
has an important impact on the unsteady aerodynamic response at a blade surface. 


Unsteady pressure-difference distributions: 
0 = 45°, G = 1.0, M_oo = 0.9, O = 180°. 
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SUBSONIC/TRANSONIC FLOW 


Finally, predictions of the aerodynamic moment versus interblade phase angle are shown for 
unit frequency torsional blade vibrations for the flat-plate and DCA cascades. Here M_ ^ = 0.9 
and those angles at which an acoustic resonance occurs are indicated by the arrows at the top of 
each figure. The unit-frequency torsional blade motions of the two cascades are stable (i.e., m/ < 0 
for all a ), but the behaviors of the flat-plate and DCA moment responses vs. interblade phase 
angle are quite different. These differences occur not only because the mean flow through the flat- 
plate cascade is entirely subsonic while that through the DCA cascade is transonic with a shock 
discontinuity, but also because of the substantial difference between the exit Mach numbers for the 
flat-plate (M 0 0 = 0.9) and the DCA (M 0 0 = 0.64) cascades. This difference implies that the two 
cascades operate in very different far-downstream acoustic response environments over almost the 
entire range of interblade phase angles. 


Unsteady moment vs. interblade phase angle; m R , m, 

0 = 45°, G = 1 , M_ oo = 0.9, cu = 1.0. 



o (deg) o (deg) 
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SUMMARY 


The linearized unsteady aerodynamic theory outlined above accounts for the effects of real blade 
geometry, mean blade loading and operation at transonic Mach numbers on the unsteady aerody- 
namic response produced by the blades of an isolated two-dimensional cascade. This theory has 
been developed to meet the requirements of turbomachinery aeroelastic designers, but it should 
also be useful for aeroacoustic design applications. The unsteady flow is regarded as a small pertur- 
bation of a fully nonuniform isentropic and irrotational mean or steady flow, which is produced by 
small-amplitude temporally and spatially (in the cascade direction) periodic structural (blade mo- 
tions) and external aerodynamic (incident entropic, vortical and acoustic disturbances) excitations. 
Thus the steady flow is determined as a solution of a full-potential boundary-value problem and the 
linearized unsteady flow as a solution of a time-independent, linear, variable-coefficient, boundary- 
value problem in which the variable coefficients depend on the underlying mean or steady flow. 

Response predictions have been presented for the blades of compressor- and fan-type cascades 
undergoing pure torsional motions. In these examples there are no incident entropy or rotational 
velocity fluctuations and, therefore, only a single field equation must be solved to determine the 
linearized unsteady flow field. The numerical results demonstrate, to some extent, the status of nu- 
merical field methods for solving the nonlinear steady and the linear, variable-coefficient, unsteady, 
boundary-value problems and illustrate partially the effects of blade geometry, mean incidence, 
shock phenomena and differences between inlet and exit free-stream conditions on the unsteady 
response at blade surfaces. 


• Linearized unsteady aerodynamic analysis 
e Effects of: 

• Blade geometry 

• Blade loading 

• Shocks and their motions 

• High frequency unsteady motions 

e Blade flutter prediction 
e Subsonic /transonic Mach numbers 
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FUTURE DIRECTIONS 


Linearizations relative to nonuniform steady flows offer great potential for meeting the needs 
of aeroelastic (or aeroacoustic) designers for efficient unsteady aerodynamic analyses that contain 
much of the essential physics associated with turbomachinery flow fields. However, before this po- 
tential can be fully realized, significant improvements in numerical solution methods for both the 
steady and linearized unsteady flows must be achieved so that reliable response information can 
be provided over the wide range of geometric configurations and flow conditions at which blade 
vibrations are of practical concern. In particular, unsteady aerodynamic analyses intended for tur- 
bomachinery aeroelastic predictions must be applicable to fan, compressor apd turbine cascades, 
low subsonic through low supersonic Mach number operation and moderate through high frequency 
structural and external aerodynamic excitations. Some needed capabilities include the ability to 
predict transonic flows (i.e., subsonic flows with imbedded supersonic regions) through fan and 
compressor cascades operating at high positive or negative mean incidence, supersonic flows with 
complicated moving shock patterns and the high frequency unsteady flows caused by incident exter- 
nal aerodynamic disturbances. Major advances in our ability to predict turbomachinery aeroelastic 
and aeroacoustic behavior should result if future research is directed toward including the effects 
of strong viscid/inviscid interactions and possibly largje-scale flow-separations within a linearized 
unsteady aerodynamic framework. Ultimately, linearized analyses which account for nonuniform 
steady flow and viscid/inviscid interaction phenomena must be extended to treat three-dimensional 
flows. 


• Subsonic/transonic flows at high incidence 


• Forced aerodynamic excitations 

•^>=f(v?s) 

• Supersonic Mach numbers 

• Viscous separation phenomena 
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